Molecular motor that never steps backwards 
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We investigate the dynamics of a classical particle in a one- dimensional two-wave potential 
composed of two periodic potentials, that are time-independent and of the same amplitude and 
periodicity. One of the periodic potentials is externally driven and performs a translational motion 
with respect to the other. It is shown that if one of the potentials is of the ratchet type, translation 
of the potential in a given direction leads to motion of the particle in the same direction, whereas 
translation in the opposite direction leaves the particle localized at its original location. Moreover, 
even if the translation is random, but still has a finite velocity, an efficient directed transport of the 
particle occurs. 
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A particle subject to a spatially asymmetric but on 
large scale homogeneous potential displays a symmetric 
diffusive motion, since the sole violation of the x — > — x 
symmetry is not sufficient to cause a net directional 
transport. As already noted more than 100 years ago by 
Curie the additional breaking of time reversal t — > — t 
symmetry [e.g. by dissipation] may lead to a macroscopic 
net velocity, so that in this case directed motion can 
result in the absence of any external force. Such sys- 
tems, known as thermal ratchets |2J, have been subject 
of much activity, both theoretical [J3| — [L 8|| and experimen- 
tal [jl9|-^5f, partly motivated by possible applicability to 
biological motors |26|-|2"s|| . 

In this Letter we study the classical dynamics of a par- 
ticle in a one-dimensional two-wave potential. The total 
potential is composed of two periodic potentials, that are 
time-independent and of equal amplitudes and periodic- 
ities. One of the potentials is externally driven perform- 
ing a translational motion with respect to the other. It 
is shown that if, in addition to the broken time reversal 
symmetry, the spatial symmetry is broken for one of the 
potentials, the relative translation can result in a two- 
fold behavior: (i) Translation in one direction causes a 
deterministic motion of the particle in the same direction, 
whereas (ii) translation in the opposite direction leaves 
the particle localized at its original location. Thus, the 
total potential acts as a ratchet in the original sense. 
Moreover, an efficient directed transport occurs even if 
the translation is random but still has a finite velocity. 
The reason for the directed transport is the existance of 
points of irreversibility in the particle trajectory. The 
high rate transport stems from the fact that if the parti- 
cle once gains a distance which is an integer multiple of 
the potential period, this distance is preserved, different 
from former ratchet systems driven by random fluctua- 
tions. 

We consider a simple ratchet type potential U(x), 
which is assumed to be continuous but not necessarily dif- 
ferentiable. It has a periodicity 6, so that H(x+b) = H(x) 
Vcc, an amplitude IIo = maxll(a;) = — minll(a;), and one 
minimum is located at x = 0, i.e. 11(0) = —IIo. We 
also assume that the potential II(x) has only one mini- 
mum and one maximum per period b, so that dH(x)/dx 



changes sign only twice for x € [0,6), cf. p9[ |. The to- 
tal potential V(x ) / y) is composed of two, not necessar- 
ily identical potentials 11(2;) and II' (a;), i.e. V(x,j) = 
Tl(x) + H'(x — 7), where 7 defines the translation. Be- 
cause of the periodicity of the potentials H(x), the po- 
tential V(x,"f) is periodic in both arguments, so that 
V(x, 7 + 6) = V(x, 7) V7. In this potential landscape, the 
deterministic equation of motion of a particle of mass m 
reads as 



r]X 



dV{x,i) 
dx 







(1) 



where the damping is denoted by rj. It should be em- 
phasized that, through the translation, energy is being 
permanently fed into the system, and that this energy 
has to be dissipated, i.e. 77 > 0. Otherwise the particle 
will gain energy until it decouples from the potential. 

First, we restrict ourselves to the case characterized by 
an overdamped motion with r)/[(2ir/b) \/mTlo} ^> 1 and 
by a slow translation with \j\/[2ny/$>Q/m] <C 1, where 7 
is the translation velocity, and relax these restrictions to- 
wards the end. In this limit, as the translation 7 is varied, 
the particle either (i) moves slowly remaining at the local 
potential minimum or (ii) if the minimum ceases to ex- 
ists, it jumps to the next minimum following the potential 
slope. The latter happens instantaneously on the time 
scale of the translation. In order to obtain the observ- 
ables, which are the trajectory x and the average velocity 
x, it is therefore sufficient to study the behavior of the 
total potential V(x, 7) and evaluate the positions of the 
minima. Below this limit is refered as 'quasistatic'. The 
initial conditions at t — are chosen as x — and 7 = 0, 
so that the particle is located at a potential minimum. 
In what follows, we use the abbreviations x = x/b mod 1 
and 7 = 7/6 mod 1. 

For simplicity, we start the discussion with a particular 
example for the potential II(x), 



n 5 (x) = n 
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which is piecewise linear, although the arguments below 
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FIG. 1. Plot of the ratchet given by eq. (Q) with asymme- 
try parameters £ = 2/5 [solid line] and f = 1/2 [dashed line, 
symmetric case]. 

apply analogously for other ratchet type potentials as 
well. The parameter £ G (0, 1) determines the asymme- 
try of the potential, with £ = 1/2 being the symmetric 
case. In Fig. [j] two realizations are shown for £ = 2/5 
and £' = 1/2. In order to introduce only one asymmetry, 
we choose the translated potential 11^, (x) to be symmet- 
ric, i.e. £' = 1/2. For the other potential IIj(x) we use 
£ < 1/2, since the case £ > 1/2 can be mapped on the 
case £ < 1/2 by replacing the asymmetry £ by 1 — £ and 
the translation 7 by —7. The case £ = £' = 1/2, where 
both potentials are symmetric, is excluded, since then the 
total potential V^i (x, 7) becomes piecewise flat, and the 
'quasistatic' treatment of eq. (0) is no longer valid [j30| . 

Let us first discuss a translation with a constant trans- 
lation velocity 7 = const, so that 7 = jt, where the 
velocity can be either 7 < or 7 > 0. In Fig. || shown 
is one cycle of a translation by —b [7 < 0] and b [7 > 0] 
for an example with £ = 2/5 and £' = 1/2. In the case 
7 < [the open circles in Fig. ||, the time evolves from (h) 
to (a)] the particle moves a distance —b, whereas in the 
case 7 > [the full circles in Fig. g, the time evolves from 
(a) to (h)] the particle, although moving locally, returns 
to its starting point. After that, the whole cycle starts 
over again. Hence, the particle moves either with average 
velocity x = 7 for 7 < or x = for 7 > 0. For the op- 
posite asymmetry £ > 1/2, the particle remains at x — 
for 7 < and moves with average velocity x — 7 > in 
the case of an opposite translation. 

This behavior can be generally understood from the 
properties of the total potential V(x, 7). In Fig. |^ shown 
are the potential V(x, 7) and, for both 7 < [left side] 
and 7 > [right side], the particle position x for a par- 
ticle located at x = for translation 7 = 0, as it fol- 
lows the changes of the potential V(x, 7). Due to the 
asymmetry of one of the constituting potentials, there 
are points of instability X in the (x, 7) plane, where a 
local minimum of the potential V(x, 7) ceases to exist. 
Hence, when the minimum disappears, a particle located 
at such a point performs an irreversible motion jumping 
to the next minimum, see Figs. p](c),(f) and Fig. ||. In 
the example considered here, this new minimum moves 
in the direction of the jump and a net transport oc- 
curs in the case 7 < 0, whereas in the case 7 > it 
moves in the opposite direction and cancels the distance 
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FIG. 2. Time evolution of the total potential V^/ (x,^f) for 
the ratchet given by eq. (^). In parallel, the respective posi- 
tions of the particle are shown, both for 7 > [full circles, 
the time evolves from (a) to (h)] and 7 < [open circles, the 
time evolves in the opposite direction from (h) to (a)]. The 
parameters are £ = 2/5 and £' = 1/2, and snapshoots are 
taken at (a) 7 = 0, (b) 7 = 1/5, (c) 7 = 2/5, (d) 7 = 41/100, 
(e) 7 = 49/100, (f) 7 = 1/2, (g) 7 = 4/5, and (h) 7 = 0. The 
arrows indicate the direction of irreversible motion of the par- 
ticle which occurs between snapshots (e) and (f) [full circle] 
and (d) and (c) [open circle]. 

gained by the jump. Under our restrictions on Tl(x), 
in particular due to the equality of the potential am- 
plitudes, all points (x, 7) € T satisfy V(x, 7) = 0. Let 
C = {(x, ; y)\V(x,j) = 0} D T be the set of all pairs (x, 7) 
for which the potential V{x, 7) = 0. Since the potentials 
II(x) are continuous, the topology of C is such that it 
consists of connected points, that form paths and inter- 
sections, see Fig. ||. The intersections correspond to the 
points of irreversibility X in which a minimum of V(x, 7) 
with respect to x ceases to exist. For the particular choice 
of the potentials H(x) given by eq. (^) one obtains the two 
points l u , = {(£,7)1(5,7) = (0, 1 - <f ), (£,£)}• Within 
the 'quasistatic' limit, jumps in the direction given by 
the potentials' asymmetry occur if the particle reaches 
these points pi] ]. 

In order to decide in general, if for a given choice of 
potentials n(x) a directed transport is possible within 
the 'quasistatic' limit, one has to examine the set C 
around the intersection points X. In Fig. [| the pos- 
sible scenarios are shown, each with sketched 'horizon- 
tal' and 'vertical' lines corresponding to V(x, 7) = 
around the intersection points. Concerning the ques- 
tion of irreversible jumps, namely if the particle reaches 
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FIG. 3. Contour plot of the total potential Vjf/(a:,7) for 
the ratchet given by eq. (^) with £ = 2/5 and £' = 1/2, the 
solid equipotential lines are placed at V^^/(x,j) = ±nIIo/5 
with 1 < n < 10 integer, and the dashed-dotted equipoten- 
tial lines indicate V( (i(x,'y) = 0. The respective trajecto- 
ries of a particle starting at position x = at translation 
7 = are shown for both 7 < [left side] and 7 > 
[right side] with thick lines, which are solid for the part of 
the trajectory where the particle remains in the minimum, 
and dashed for the irreversible jumps. The arrows indi- 
cate the time development for 7 < [downward arrows] 
and 7 > [upward arrows]. The points of irreversibil- 
ity T = {(55,7)102,7) = (0,1/2), (2/5, 2/5)} are marked by 
dashed circles. 

the points of irreversibility as the translation is 
monotonously varied, one has to examine the behavior 
of the horizontal line, i.e. how this line is bent at the 
intersection point with respect to the direction of the 
translation. In the upper left part of Fig. || the hori- 
zontal line is bent downwards towards smaller values of 
7 and hence opposite to the direction of the translation 
7 > on both sides of the intersection point. This means 
that (i) the minimum in which the sketched particle is 
located moves towards the intersection point as 7 is in- 
creased [because of the downward bending on the parti- 
cle's side], and (ii) the minimum ceases to exist at the 
intersection point with a local slope such that the par- 
ticle performs a jump leftwards to smaller values of x 
[because of the downward bending on the side opposite 
to the paricle]. In the lower left part of Fig. [| the hor- 
izontal line is bent upwards on the particle's side of the 
intersection point, so that the minimum does not move 
towards the intersection as 7 is increased, but remains 
always right to it at larger values of x. This behavior 
is independent of the bending on the side of the inter- 
section point opposite to the particle, either upward, as 
shown in the figure, or downward. In the upper right 
part of Fig. ^ a third theoretically possible topology is 
shown. In this case, the particle reaches the intersection 
point [because of the downward bending on the particle's 
side] , but the intersection point is not a point where the 
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FIG. 4. Possible topologies for the set C around the points 
of irreversibility I. The dashed lines indicate V{x, 7) = 0, and 
the areas with V(x, 7) > and V(x, 7) < are indicated. The 
particle's position is shown as full circle. 

minimum ceases to exist [because of the upward bend- 
ing on the side opposite to the particle]. However, this 
topology cannot occur under our restrictions on the po- 
tentials Tl(x). For the potentials H(x) given by eq. (^) 
and for £ 7^ 1/2 or £' 7^ 1/2 the resulting topology of C is 
always such that the points of irreversibility are reached 
within the 'quasistatic' limit. Hence, depending on the 
asymmetry of the potentials, one observes transport for 
either 7 < or 7 > 0, and no transport in the case of an 
opposite translation. 

We can extend the model beyond the translation with a 
constant velocity, by assuming for instance another sim- 
ple scenario with an oscillatory translation of the form 
7 = 70 + 71 sin(27ra->t) with a driving frequency u> and 
77(71 |w&/no <C 1. In this case, one finds an average ve- 
locity x = ±nbuj with n > integer, where the sign 
depends on the asymmetry of the potentials IL(x). The 
actual value of n depends on the offset 70 and on the am- 
plitude 71 , since these values determine how many times 
the points of irreversibility I are reached during one cy- 
cle. For the potentials H(x) given by eq. ^ with £ = 2/5 
and £ = 1/2, a choice of 7o = 456/100 and 71 = 6/20 
results in an average velocity of x = ~buj. This means 
that in order to make the particle gain a distance 6, the 
translation has only to be varied twice by 6/10. 

As a third possible scenario for the translation 7, we 
assume 7 to vary randomly by following the trace of a 
random walker y that has locally a finite constant velocity 
< \y\ = const, but an average velocity y = 0. For this 
scenario we relax the restriction of overdamped motion to 
demonstrate the validity and accuracy of the 'quasistatic' 
treatment, and numerically integrate eq. (Q) with a finite 
damping and a finite translation velocity. However, one 
has to keep the restriction |7|/[27r^/ < I>o/m] <C 1, since 
otherwise the changes in the potential due to the trans- 
lation are too fast for the particle to follow, and the par- 
ticle decouples from the potential. In Fig. || a realization 
of a random walker trace y and the resulting particle tra- 
jectory x for £ = 2/5 and £' = 1/2 are shown. We would 
like to emphasize that, despite the relatively small damp- 
ing, the numerically obtained trajectory differs from the 
one obtained under overdamped conditions only by small 
oscillations. Although the random walker, and hence the 



3 




??/[(27r/6)Vmi o] = 0. 2. 
locity |y|/(27r A /*o/m) 



FIG. 5. Plot of a random walk trace y [(a), in grey] 
used as translation and the numerically obtained parti- 
cle trajectory x [(b), in black] for asymmetry parame- 
ters £ = 2/5 and {' = 1/2, and dissipation constant 
The random walker moves with ve- 
= 0.02 and chooses the direction anew 

every t/r = 1 with r = [(27r/&)y^&o/m] _1 , thus gaining a 
distance \Ay\/b = 0.02 during each step. 

translation, have an average velocity y — resp. 7 = 0, 
the particle's average velocity x is nonzero. The sign of 
the velocity depends on the asymmetry of the potentials 
II(a;), for our choice of asymmetry one finds x < as 
expected. The transport is much more efficient than in 
former ratchet systems driven by random fluctuations. 
Here, if the particle gains a distance d(t) = \x(t)\ equal 
to an integer multiple of the potential period b at a time 
tm d(t n ) = nb, the distance will never become smaller 
again, and d(t) > nb^t > t n . Putting this in the context 
of molecular motors, it means that the suggested molec- 
ular motor never executes a step backwards. These back- 
ward steps usually limit the efficiency of the motor ||. 

Different realizations of the model can be thought of. 
One possibility is to put a small particle in a potential 
created by an optical tweezer such as in pp| , but with 
two time-independent potentials added on top of each 
other at the same place and having a certain phase shift. 
Although the phase shift is changed randomly, but with 
a finite velocity, an efficient transport of the particle is 
predicted in a certain direction determined by the asym- 
metry of one of the potentials. It is possible, for instance, 
to induce the translation of the potential through the 
coupling of a varying internal degree of freedom of the 
particle to the non-translated potential 0] . This opens 
various possibilities for the construction of micro- and 
nanoscale devices such as pumps and motors based on 
the presented ratchet type system. 
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